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Abstract 

We propose a systematic procedure for constructing effective models of strongly correlated ma- 
terials. The parameters, in particular the on-site screened Coulomb interaction U, are calculated 
from first principles, using the GW approximation. We derive an expression for the frequency- 
dependent U(u) and show that its high frequency part has significant influence on the spectral 
functions. We propose a scheme for taking into account the energy dependence of U(cj), so that a 
model with an energy- independent local interaction can still be used for low-energy properties. 



1 



I. INTRODUCTION 



Lattice fermion models such as the Hubbard model or the Anderson impurity model and 
their extensions have played a major role in studying electron correlations in systems with 
strong onsite correlations. Despite the widespread use of these models, little justification has 
been given in using them. The models are postulated on the basis of physical intuition. In 
particular, the models employ parameters, such as the famous Hubbard interaction U, which 
are normally adjusted to serve the given problem. Without judicious choice of parameters, 
the model may yield misleading results or in the worst case, the model itself is not sufficient 
to describe the real system. One can define rigorously these concepts in the path integral 
formulation of the many body problem by performing a partial trace over the degrees of 
freedom that one wants to eliminate, and ignoring the retardation in the interactions gen- 
erated by this procedure. However, this elimination of the degrees of freedom is very hard 
to perform for real materials. It is therefore very desirable to figure out a systematic way of 
constructing low energy effective models with well defined parameters calculated from first 
principles such that the model can quantitatively reproduce and predict physical proper- 
ties of interest of the corresponding real system, especially when the correlation effects are 
crucial. 

Another important issue that has not received sufficient attention is the role of energy 
dependence of the screened local Coulomb interaction U. Model studies investigating the 
lmp o rtm ce „ thig h-e Mrgy states in the H ubbard mod e, c m b e foued in flflfl. A dynamic 
Hubbard model has also been considered. |4|. In most cases however U is assumed to be 
static, but on the other hand we know that at high energy the screening becomes weaker 
and eventually the interaction approaches the large bare Coulomb value, which is an order 
of magnitude larger than the static screened value. Of course the high-energy part of the 
Coulomb interaction has in some way been down-folded into the Hubbard U but it is not 
clear how this downfolding is actually accomplished. 

A number of authors have addressed the problem of determining the Hubbard U from 
first principles. One of the earliest works is the constrained local density approximation 
(LDA) approach 0] where the Hubbard U is calculated from the total energy variation 
with respect to the occupation number of the localized orbitals. An approach based on the 
random-phase approximation (RPA) was later introduced p, which allows for the calcula- 
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tions of the matrix elements of the Hubbard U and its energy dependence. This was followed 
by a more refined approach for calculating U. Q. A yet different approach, computes the 
matrix elements of the Coulomb interactions screened in real space and assumes a Yukawa 
form to extract the Hubbard U and the other interactions which determine the multiplet 

n 

splittings. [8fl. 

The purpose of the present work is to develop a precise formulation for a systematic con- 
struction of effective models where the parameters are obtained from realistic first-principles 
electronic structure calculations. In particular, we concentrate on the calculation of the Hub- 
bard U and demonstrate the importance of its energy dependence. We show that a static 
Hubbard hamiltonian, obtained from a naive construction in which this energy dependence 
is simply neglected, fails even at low energy. This static model can be appropriately mod- 
ified however, by taking into account the feedback of the high-energy part of U into the 
one-particle propagator. We illustrate our scheme in transition metals, concentrating on Ni 
as an example, since it is a prototype system consisting of a narrow 3d band embedded in a 
wide band. Furthermore, Ni is one of the most problematic case from the viewpoint of the 
LDA. 

II. THEORY 

Let us suppose that the bandstructure of a given solid can be separated into a narrow 
band near the Fermi level and the rest, like, for example, in transition metals or 4f metals. 
Our aim is to construct an effective model which only includes the narrow 3d or 4f band. 
The effective interaction between the 3d electrons in the Hubbard model can be formally 
constructed as follows. We first divide the complete Hilbert space into the Hubbard space 
{ipd}, consisting of the 3d states or the localized states, and the rest. The bare Green's 
function G d , spanning the d-subspace is given by: 
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Let P be the total (bare) polarization, including the transitions between all bands: 

occ unocc 

P(v,v';u) = J2H ^(r)^(r')^(r)^(r') 

i j 

1 1 

a; — £j + Si + i0+ u; + Sj — Si — i0 + 

P can be divided into: P = P d + P r , in which includes only 3d to 3d transitions (i.e 
limiting the summations in (J2J) to i,j G {V'd})) an d P r be the rest of the polarization. The 
screened interaction W on the RPA level is given by 
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where we have defined a screened interaction W r that does not include the polarization from 
the 3d-3d transitions: 

W r (u) = [l-vP r {u)]- l v (4) 

(we have not explicitly indicated spatial coordinates in this equation). The identity in 
(JHJ) explicitly shows that the interaction between the 3d electrons is given by a frequency- 
dependent interaction W r . It fits well with the usual physical argument that the remaining 
screening channels in the Hubbard model associated with the 3d electrons, represented by 
the 3d-3d polarization P d , further screen W r to give the fully screened interaction W. 

We now choose a basis of Wannier functions {<f)R n }, centered about atomic positions R, 
corresponding to the 3d Bloch functions {V'kn}, and consider the matrix elements of the 
(partially screened) frequency- dependent Coulomb interaction W r : 

U Rin R 2 n',R 3 mR 4 m'(r -r') = J d 3 rc?V^ iri (r)^^ (r) W r {r, r';r - r , )(j) R3m {r') (j)^, (r') (5) 

We would like to obtain an effective model for the 3d degrees of freedom. Because of the 
frequency dependence of the U's (corresponding to a retarded interaction), this effective 
theory will not take a hamiltonian form. We can however, write such a representation in 
the functional integral formalism^] by considering the effective action for the 3d degrees of 



freedom given by: 

S = J drdr' [- £ dUr) R i R/r Ar ~ r')d Rnl {r') 

: d R in (r)d R2n/ (T) : U RinR2n > jRamRim ,(T - r') : d R3m {T r )d R4m >(T') : 



(6) 



where : d^d : denotes normal ordering, which accounts for the Hartree term, and the sum- 
mation is over repeated indices. When using a Wannier transformation which does not mix 
the d-subspace with other bands, the Green's function can be taken, to first approximation, 
to be the bare Green's function G° dd constructed from the Bloch eigenvalues and eigenfunc- 
tions. If instead an LMTO formalism |l0| is used, one should in principle obtain from a 
downfolding procedure onto the d-subspace, i.e perform a partial trace over s,p degrees of 
freedom (e.g to first order: _1 = [G^] -1 - [G ds ]' l G ss [G sd ]' 1 ). 

In the following, we retain only the local components of the effective interaction on the 
same atomic site. This is expected to be a reasonable approximation because the 3d states 
are rather localized. The formalism may be easily extended to include intersite Coulomb 
interactions if necessary. Hence, we consider the frequency-dependent Hubbard interactions: 

U nn r, mm ,{r -t') = J rf 3 rrf 3 rV;(r)0„,(r)H/ r (r, r'; r - r')C(r')<Mr') (7) 

with (p n being the Wannier orbital for R=0. In order to illustrate the procedure within 
the linear muffin-tin orbital (LMTO) basis set, we use instead of the Wannier orbital the 
normalized function head of the LMTO <pH which is a solution to the radial Shrodinger 
equation matching to a Hankel function at zero energy at the atomic sphere boundary. 

In this paper, we investigate the importance of the energy- dependence of U. Therefore, 
we shall compare the results obtained from (jUJ) with those of a Hamiltonian approach in 
which one would construct a Hubbard model with a static interaction U : 

H = c Rn h Rri:R i n ic R i n i + - c Rn c Rn 'U nn > , mm >c Rm c Rrn ' (8) 

Rn,R'n' R,nn',mm' 

It seems natural to identify the static Hubbard U with the (partially) screened local inter- 
action in the low- frequency limit W r (u = 0). Note that the Hubbard model (jHJ) has been 
constructed in the most naive manner, by simply taking the quadratic part to be the d-block 
of the non-interacting hamiltonian. 
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FIG. 1: respect to the \$E\_g\$ channel. Other channels give almost the same result. 

In order to compare the results obtained from the full dynamical U to those of the 
static Hubbard model, we need to solve (jBJ) and (jBJ) within some consistent approximation 
scheme. In the following, we adopt a "GW universe" where the exact self-energy for the 
solid is assumed to be given by the GW approximation (GWA) [lllQl, and where the GW 
approximation is also assumed to be a reliable tool in solving the effective models © and 
(JHJ). This allows us to make a proper comparison between the "exact" self-energy and the 
Hubbard model self-energy. If the assumption of static U is valid, the Hubbard self-energy 
and the true self-energy (both within the GWA) should be close to each other, at least for 
small energies. Or equivalently, the spectral function for small energies should resemble that 
of the full one. 

III. RESULTS AND DISCUSSIONS 
A. Comparing self-energies 

The screened interaction with and without the 3d-3d transitions is shown in Fig. Q in 
the case of Nickel. Here and in all the following, a spin-unpolarized (paramagnetic) solution 
is considered. At low energies, the (partially) screened interaction W r without the 3d-3d 
transitions is larger than the full one W, and at high energies they approach each other, as 
anticipated. Related calculations have also been performed by Kotani 

We first compare the self-energy obtained from a GW treatment of the full system, given 
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FIG. 2: The self-energy of the real system (nickel) and the Hubbard model with a frequency- 
dependent U. 



to the self-energy obtained from the effective model © with an energy- dependent interaction 
U{u) = W r {u). Because of (jHJ), the screened interaction corresponding to this U (u) is simply 
W, and the corresponding self-energy reads: 



The difference between this expression and the GWA for the full system (EqEJ) is that in 
(IK)jl only the d-block of the Green's function has been included (since the effective action 
was written for the d-band only). Hence, the two self-energies differ by a term G r W, with 
G r = G — Gd- We expect that the wavefunction overlap between two 3d states (one from G<i 
and the other from the 3d state appearing in the matrix element of S<j) and other non-3d 
states is small so that should be close to the true E. In Fig. |2l the two self-energies are 
displayed (more precisely, in this figure and in all the following, we display the matrix element 
of the self-energy in the lowest 3d state (band number 2), at the T-point, corresponding to 
an LDA eigenenergy -1.79 eV). We observe that the two self-energies indeed almost coincide 
with each other, even at high energies. Hence, we conclude that the effective Hubbard model 
(JHJ) for the d-subspace, with an energy-dependent interaction, provides a reliable description 
of the real system. 
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FIG. 3: The real parts of the self-energy of the real system (solid) and the Hubbard model (dash) 
in the GWA. 

We now turn to the self-energy associated with the static Hubbard model (JEJ) with U = 
W r (0) ~ 3.5eV. The local screened interaction for this model, within the GW approximation, 
is given by: 

W d {u) = [l-UP d (uj)}- l U (11) 
and the self-energy of the Hubbard model in the GWA thus reads: 

Ef (r, r'; u) = ^ / dJe^G A {r, r'; u + u/)W d (r, r'; u') (12) 
Ztt J 

Note that the difference between this static Hubbard model self-energy and that of the 
effective model with a frequency-dependent interaction (JHJ) relies in the use of a different 
form of the screened interaction W d instead of the full W. In Fig. (|5|L the real part of this 
self-energy with that of the full GWA self-energy for nickel are shown. Since the energy 
scale of the self-energy of the real system is determined by the bare Coulomb interaction v 
whereas the Hubbard self-energy is set by U, the latter has been shifted so that it is equal to 
the former at the LDA eigenvalue (-1.79 eV) of the band we have considered, at the T-point. 
The difference in magnitude of the self-energies is not important since it simply shifts the 
spectrum (or, said differently, we have compared differences E — /i from the values of the 
chemical potential obtained in the various schemes). However, the difference in the variation 
of the self-energy with respect to energy matters, since it will give a different quasiparticle 
weight Z — [1 — dKe'E/duj]~ 1 and affect the spectral function. As can be deduced from 
the figure, the Z factor of the Hubbard model taken at the energy of the quasiparticle band 
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FIG. 4: The imaginary parts of the self-energy of the real systems (solid) and the Hubbard model 
(dash) in the GWA. 

at the T-point (~ — 1.79eV) is much closer to unity as compared to the true (full GW) 
one because the former already contains the renormalization from the plasmon. Hence, 
neglecting frequency- dependence directly affects the physical results, even in the low- energy 
range. We shall see below, however, that it is possible to modify the static model in such a 
way that an accurate approximation is obtained at low energy. 

In Fig. |3] the imaginary part of the self-energy is shown. Here we see that Im T>f for 
the Hubbard model is peaked around 5 eV since there are no states above or below the 
3d band. As a consequence, the real part of the self-energy exhibits the Kramers-Kronig 
behavior at around -5 eV. Within an energy region spanning about twice the 3d bandwidth, 
the imaginary part of the Hubbard model self-energy, in contrast to the real part, is not in 
bad agreement with the full one. 

These findings can be understood qualitatively by considering more explicit expressions 
of the self-energy obtained within the GW approximation, for an effective model of the 
d-subspace defined on the periodic lattice. The imaginary part reads: 

ImE^u,) = ^ImW/ nm (g,^ - e^ q ) [n F (e%_ q ) + n B {u - e£.J (13) 

q,m 

In this expression, n, m are band indices, ejj corresponds to the n-th non-interacting band and 
np (resp. rig) is the Fermi (resp. Bose) function. ImW is the spectral function associated 
with the effective interaction (to be taken as W if the self-energy of the frequency- dependent 
effective model is considered, and as Wd if that of the static Hubbard model is considered). 
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FIG. 5: The spectral functions of the real system (solid) and the Hubbard model (dash) in the 
GWA. 

From this expression, one sees that if one considers an energy u, the bands contained in the 
energy interval [0, \u\ ] are the only ones contributing significantly to ImS. This is the reason 
why the imaginary part of the self-energy at low-energy will be correctly reproduced by the 
effective low-energy model (provided of course that the spectral function ImW is correctly 
approximated at low energy). In contrast, the real part of the self-energy is obtained from 
the Kramers- Kronig relation in the form: 

ReE n (k,u) = - l -P f + °°du ^WUa, ^^-^^ (14) 

J 00 q,m K-q 

Because the principal- part integral extends over the whole frequency range, high-frequency 
contributions influence the self-energy even at low frequency. As a result, an accurate de- 
scription of the real part of the self-energy cannot be obtained within the naive construction 
of the static Hubbard model because the effective interaction is not correctly approximated 
over the whole frequency range. This formula also suggests a way to appropriately mod- 
ify the effective static model in order to obtain an accurate description at low energy, as 
discussed below. 

B. The puzzle of the satellite 

Fig. El compares the spectral function obtained for the full system within the GWA, and 
that of the static Hubbard model (jHJ). A striking difference between these two results is the 
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absence of the 6 eV satellite in the GWA for the full system, while the static Hubbard model 
displays a satellite feature. This is to be expected from the structure of the self-energy. The 
GWA self-energy of the full system continues to grow and reaches a maximum around the 
plasmon excitation at around 25-30 eV, while the self-energy of the static Hubbard model 
has a maximum around the width of the 3d band (4 eV). In comparison with the true 
plasmon, the "plasmon" of the Hubbard model has a much lower energy, of the order of 
U, which results in the " 6 eV satellite" . One could blame the appearance of a satellite in 
our static Hubbard model calculation on the fact that we have used a (non self-consistent) 
GW approximation. However, more accurate treatments of the static Hubbard model (such 
as DMFT) do preserve this feature, which has in this context a natural interpretation as 
a lower Hubbard band. On the other hand, it could be that the experimentally observed 
Nickel satellite has a somewhat different physical origin and that the satellite obtained in 
static Hubbard model calculations is spurious in the context of Nickel (e.g because it is not 
legitimate to use a low-energy effective model in this energy range). In our view, this issue 
is an open problem which deserves further work. 

C. Improving the effective static model 

The preceding discussion shows that an effective model for the d-band can be constructed, 
which accurately reproduces the full results over an extended energy range, provided the 
energy- dependence of U (cu) is retained. However, performing calculations with an energy- 
dependent Hubbard interaction is exceedingly difficult. A more modest goal is to obtain 
an effective model which would apply to some low-energy range only (say, < A, with 
A a cutoff of the order of the d-bandwith). In order to achieve this goal, we propose to 
adopt a renormalization group point of view, in which high energies are integrated out in 
a systematic way. Following this procedure, an appropriate low-energy model with a static 
U can be appropriately constructed. As we shall see, the bare Green's function defining 
this low-energy model does not coincide with the non-interacting Green's function in the 
d-subspace (we have seen that this does not lead to a satisfactory description, even at low 
energy) . 

Let us illustrate this idea within the GW approximation. The full Hilbert space is divided 
into the Hubbard space, comprising the 3d orbitals, and the downfolded space, comprising 
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the rest of the Hilbert space. This approach is complementary to the one in , where 

the division is done in real space (onsite and off sites). Within the GWA we may write the 
full self-energy as follows: 

GW = G d W d + G d (W - W d ) + G r W. (15) 

where G = G d + G r . The term G d (W — W d ) represents the high-energy contribution of 
the screened interaction. This is the main source of error in the naive static limit, as 
discussed in the previous sections. The term G r W is not obtainable within the Hubbard 
model, even when a frequency-dependent U is employed since G r resides in the downfolded 
space. This term was shown to be small. Its effect at low-energy can also be taken into 
account by appropriate modifications of the one-particle propagator, but we shall neglect it 
for simplicity. 

We consider the following Hubbard model with a static U = W r (0), but a modified 
one-particle propagator Go, defined by the action: 

Sh = -J drdr'J2dUr)[Gl n , m A-\r-ry Rn/ (r') 

+ \J dT ^l'- d Rm( T ) d R2n'(r) : U RinR2n ^ RsmRim , : dJ feTO (r)dfl 4m /(r) : (16) 
The self-energy of this static Hubbard model in the GWA is: 

Xf = G d W d (17) 

where the new effective interaction is W d = U [1 — UPd] , with P d constructed from the new 
Green's functions G d (schematically P d = G d - G d ). We request that the interacting Green's 
function of this modified static model coincides with that of the Green's function calculated 
with the frequency- dependent interaction in the low energy range \uo\ < A:), that is: 

G d 1 - Ef ~ G d l - G d W d - G d (W - W d ) for \u\ < A (18) 

Using the identity G^ 1 — GW = [1 — UP]^ 1 G^ 1 , this can be rewritten as: 

[1 - UP,]- 1 Gf ~ [1 - UP,]- 1 Gf - G d (W - W d ) for M < A (19) 

This is an integral equation which determines in principle the modified bare Green's function 
G d to be used in the "downfolded" static action ()16j) . To first approximation, one can neglect 
the polarization terms in this equation, and obtain the first order modification of G d as: 

G? = G; 1 - G d {W — W d ) H (20) 
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The first correction appearing in this equation is precisely the contribution coming from the 
high-energy part of the screened interaction. We have explained above that this correction is 
not small, which is the reason for the failure of the "naive" static Hubbard model using the 
non-interacting G d - Dividing the screened interaction into a low-energy part for \u\ < A and 
a high-energy part for |u;| > A, we can use the explicit forms (|13I14|) given in the previous 
section to obtain this first correction in the form: 

1 f x-^ n F (ef J + n B (u) 

J\v\>l\ qm k-q 

In particular, we see by expanding this expression to first order in lu, that the low fre- 
quency expansion of the modified one-particle propagator reads: G d 1 (k,u + i0 + ) n = 
(1 + a/c in )co> — e£ + ■ • • . The coefficient ak, n is a partial contribution to the quasi-particle 
residue a fcj „ = -~ J2 m I\u\>a dvlmW nm (q, v) ^^"^"f — ■ In practice this integral can be 
easily evaluated using the well known plasmon pole approximation, which should contain 
most of the high energy contribution. This correction insures that the quasi- particle residue 
is obtained correctly from this "downfolded" static model. We do not, however, expect that 
this improvement solves the problem with the satellite discussed in the previous subsection. 



D. Beyond the GW approximation 

In strongly correlated systems, it is known that the GWA is not sufficient and improve- 
ment beyond the GWA is needed. Much of the short coming of the GWA probably orig- 
inates from the improper treatment of short-range correlations within the random-phase 
approximation (RPA). Thus, one would like in the first instance attempt to improve these 
short-range correlations, which are essentially captured by the Hubbard model. The formula 
()15j) suggests a natural way to do this, as follows. The contribution from the frequency- 
dependent U as well as the contribution from the downfolded Hilbert space are treated 
within the GWA. The self-energy G d W d corresponding to the GW self-energy of the Hub- 
bard model with a static U can then be replaced by that obtained from more accurate 
theories such as dynamical mean-field theory (DMFT) Q| or from exact methods such as 
Lanczos diagonalization. Thus, if we use the LDA to construct Gd, the correction to the 
LDA exchange-correlation potential reads: 

AS = Y, H + G d (W - W d ) + G r W - v xc , (22) 

13 



where E# is the Hubbard model self-energy obtained from more accurate methods, replacing 
GdWd, We note also that the scheme in (fToT) avoids the problem of double counting, inherent 

n n 

in LDA+U [15] or LDA+DMFT methods. [19] In this way, it is possible to calculate the 
Hubbard U using the response function constructed from the LDA bandstructure. 

The application of DMFT requires the Hubbard U for the impurity model. To calculate 
the Hubbard U for an impurity model, it is necessary to downfold contributions to the po- 
larization P from the neighboring sites. Here, the identity in Eq. shows its usefullness. 
Since the formulation in (j3J) is quite general, we merely need to redefine Pd to be the onsite 
polarizations or equivalently, we redefine P r to also include polarizations from the neigh- 
boring sites. The computational result for Ni shows that the difference between the lattice 
and the impurity Hubbard U is rather small, essentially negligible. Somewhat different but 
related calculations of the Hubbard U for the impurity model has also been performed in 
[3] for Fe and Ni, which confirm the result obtained in our formulation. 

Another feasible approach for calculating physical quantities of the derived Hubbard 
model is by using the path integral renormalization group (PIRG) method. 17] One advan- 
tage of this method is the possibility of using the lattice Hubbard model as opposed to the 
impurity model, thus including spatial fluctuations and possible symmetry breaking. The 
method is also suited for studying single-particle Green's functions as well as thermody- 
namic quantities, in particular for accurate determination of the phase diagram, which often 
requires a determination of the possible symmetry breaking of the Hubbard Hamiltonian 
after taking account of spatial and temporal fluctuations on an equal footing. In calculating 
the self-energy, one may substitute G^Wd by the self-energy obtained within the correlator 
projection method Q] together with the PIRG. These form a future challenge in the field 
of strongly correlated materials. 



IV. SUMMARY AND CONCLUSION 



In conclusion, we have investigated the construction of effective models for the correlated 
orbitals in materials with a natural separation of bands. We have shown that, if one retains 
the full frequency dependence of the local components of the screened interaction, an accu- 
rate effective model can be obtained over an extended energy range. Simply neglecting this 
energy dependence and using the non-interacting hamiltonian into a static Hubbard model 
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does not provide an accurate description even at low energy. However, a proper modification 
of the bare propagator, obtained by integrating out high energies, allows for the construc- 
tion of an effective Hubbard model which describes the low-energy physics in a satisfactory 
manner. 
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